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In zero magnetic field the superconductor Sr2Ru04 is believed to have a chiral spin triplet pairing 
state in which the gap function d-vector is aligned along the crystal c-axis. Using a phenomenological 
but orbital specific description of the spin dependent electron-electron attraction and a realistic 
quantitative account of the electronic structure in the normal state we analyze the orientation of 
the spin triplet Cooper pair d-vector in response to an external c-axis magnetic field. We show 
that for suitable values of the model parameters a c-axis field of only 20 mT is able to cause a 
reorientation phase transition of the d-vector from along c to the a — b plane, in agreement with 
recent experiments. 

PACS numbers: PACS numbers: 74.70.Pq, 74.20.Rp, 74.25.Bt 



INTRODUCTION 



susceptibility has the following tensor form 



The 1.5 K superconductor S^RuCU 1] is widely be- 
lieved to be a rare example of a spin-triplet Cooper paired 
system 0, Q . Unlike heavy fermion materials which are 
possible candidates for spin triplet pairing (such as UPta) 
there are no 4/ or 5/ heavy elements in the unit cell, 
and so the effects of spin-orbit coupling should be rela- 
tively weak. In this respect paring in Si~2Ru04 should 
be closely analogous to the case of superfluid 3 He. One 
should expect that weak or moderate £>-fields would be 
able to rotate the d-vector order parameter and induce 
domain walls, textures and other topological defects. As 
in the case of 3 He, the experimental observation of such 
d-vector rotations would be a decisive test of the pairing 
symmetry. Systematic study of such rotations as a func- 
tion of various physical parameters would both confirm 
the pairing symmetry and also place strong constraints 
on the pairing mechanism and microscopic Hamiltonian 
parameters. 

One of the strongest pieces of evidence for spin triplet 
Cooper pairing in Sr2Ru04 was the observation that 
the electronic spin susceptibility Xs(T) remained con- 
stant below T c for magnetic fields in the a-b plane 0, [B| . 
This observation is inconsistent with a spin-singlet pair- 
ing state (s or d-wave) for which Xs(T) = x-nY(T), where 
Xn is the normal state Pauli spin susceptibility and Y(T) 
is the Yoshida function. On the other hand the obser- 
vations would be immediately consistent with a chiral 
symmetry paring state of the form 



d(k) = (sin fee +isinfcj / )e z , 



(1) 



1 
Xs(T)= Xn ( 1 
Y(T) 



(2) 



which below we shall refer to as pairing state (a). For 
such a pairing state is is well known [6| that the spin 



and so the susceptibility tensor is constant as a func- 
tion of T for field directions perpendicular to the d vec- 
tor. The 17 O Knight shift experiments performed in (ab) 
plane magnetic fields [J] and neutron scattering experi- 
ments [5] both observed a constant susceptibility below 
T c , a result which uniquely points to a chiral triplet par- 
ing state with d || e z . 

More recently Murakawa et al. 7] have measured 
the 101 Ru-Knight shift of Sr2Ru04 in a superconduct- 
ing state under the influence of magnetic field parallel 
to the c-axis. In contradiction to expectations from Eqs. 
(fTj) and @ they found that its value is also unchanged 
from the normal state value below T c . They remarked 
that this result would be consistent with Eq. ([2]) only if 
it is assumed that the d vector rotates away from the e z 
direction under the influence of the external field. Such a 
rotation of the d vector is well known in the case of super- 
fluid 3 He-A Q, where to minimize free energy d orients 
itself perpendicular to the external field unless pinned by 
surface effects. However, in Sr2Ru04 it is expected that 
spin-orbit coupling would fix the chiral state d vector to 
the crystal c-axis [g, 0, [lij II | • Murakawa et al. Q were 
able to measure Xs (T) down to fields as low as 20 mT at 
80 mK, showing that it was equal to Xn to within experi- 
mental accuracy at all temperatures and fields measured 
below B C 2 (75 mT for c-axis fields). The implication of 
this result, that d vector rotation must occur at below 
20 mT, therefore places very strong constraints on the 
strength of pinning by the spin-orbit coupling. 

In this work we examine the combined effects of spin- 
orbit coupling and c-axis magnetic field on the chiral 
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state of Si'2Ru04. We show that a realistic physical 
model, which we have previously shown to be consis- 
tent with a wide range of experimental data, allows us to 
find reasonable parameter ranges where the chiral state 
Eq. JT]) is stable in zero field, but where a transition to 
another state occurs even for fields of order 20 mT or 
smaller. The minimum energy states in finite c-axis field 
are not simply the chiral state with d vector rotated to 
the a — b plane, but rather are non-chiral pairing states 
of the form 



1.0 



d(k) = (sin k x , sin k y , 0) 
d(k) = (sin k y , — smk x , 0) 



(3) 



on the 7 Fermi surface sheet (which below we refer to 
as (&) and (c) respectively). The spin susceptibility for 
either of these states is of the form \v3^ 



Xs (-0 2 ^- T - 



Y{T) 
1 + Y(T) 
2 



(4) 



corresponding to constant spin susceptibility for c-axis 
fields consistent with the results of Murakawa et al. • 
Since these states are symmetry distinct from the chi- 
ral state Eq. ((T|) the transition from one to another should 
occur at a finite external field, B t , and it should be ac- 
companied by a finite entropy change AS. Below we 
map out a generic phase diagram for the transitions be- 
tween pairing states of Sr 2 Ru04 in a c-axis field, assum- 
ing parameter values consistent with the Murakawa's et 
al. experiments Q • We also estimate the entropy change 
associated with the transition, and comment on its ex- 
perimental observability. Our preliminary results have 
been presented in (l3| . 



THE MODEL 

As the pairing mechanism operating in the strontium 
ruthenate superconductor is not known 0, [3} it is im- 
portant to understand the experiments on the basis of 
semi-phenomenological models [3 ITlj . Specifically we 
write the effective pairing Hamiltonian as 



l O 2^ 2^ £ a <y<y' c ir, 



i ma - c im'a' 



i,aa k mm 



2 ^rnmJ i^j)^ima^jm' ff^jm' ^CimS i (5) 

ijmm' aft'fS 

where m and m! refer to the three Ru ti g orbitals a = d xz , 
b = dy Z and c = d xy , and i and j label the sites of a 
body centered tetragonal lattice. The hopping integrals 
tmm'(ij) and site energies e m were fitted to reproduce 
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FIG. 1: (Color online) Temperature dependence of the nor- 
malized spin susceptibility in the superconducting state for a 
number of order parameters. The curves marked (a) and (b) 
show the c-axis susceptibility for Sr2Ru04 and chiral symme- 
try state, (o), defined by Eq. [S]and Eq. [7] and the ab-plane 
susceptibility of the non-chiral triplet state, of symmetry (b), 
corresponding to Eq. [3] Curves labelled (s) and (d) show the 
Yoshida function for single band superconductor with two di- 
mensional tight binding spectrum and order parameter of s 
and d wave symmetry, respectively. 



the experimentally determined Fermi surface 16], 17 1. A 
is the effective Ru 4d spin-orbit coupling parameter, and 
the effective Hubbard parameters U^J S (ij) are gener- 
ally both spin [3, |T|| and orbital dependent [3, EE, [2(| • 

Our approach is based on self consistent solution of the 
Bogolubov-deGennes equations for the model ([5]) in the 
spin triplet channel for each possible symmetry distinct 
order parameter. We have shown elsewhere 2o|, 21 1 that a 
minimal realistic model requires two Hubbard U param- 
eters: Uu for nearest neighbor in-plane interactions be- 
tween Ru d xy orbitals (corresponding to the 7 Fermi sur- 
face sheet) and U± for out of plane nearest neighbor inter- 
actions between Ru d xz and d yz orbitals (corresponding 
to the a and (3 Fermi surface sheets). When these two 
parameters are chosen to be spin independent constants 
(for A = 0) and to give the experimental T c = 1.5K on 
all three a, (3 and 7 Fermi surface sheets, then the calcu- 
lated specific heat, thermal conductivity and a — b plane 
penetration depth are in very good agreement with ex- 
periments in mm 1 . 



Figure ([T]) shows the susceptibilities calculated within 
the model for the (a) and (b) states with the symmetries 
given by Eqs.[T]and[3] For the chiral state (a) the triplet 
d-vector is defined in all 3 bands as d(k) = A mm '(k)e z , 
with A mm <(k) denoting contributions from different or- 
bitals which are given by 



A cc (k) = A cc (T)(sin k x + i sin k y ) 



(6) 
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for the Ru c(= d xy ) orbitals and, 



,(k) 



— cos -r 
2 2 
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a, 6 corresponding to the Ru d xz and d 



interaction 
than [/ = 



for m, m' 

orbitals [2fJ]. The susceptibility results, for the chiral (a) 
state show that the dependence of \ f° r c-axis B-fields is 
is relatively structureless. The relevant Yoshida function 
Y(T) is equal to unity at T = T c and drops essentially 
linearly to zero at T = 0, consistent with the expecta- 
tion from EqJ2] and the existence of the horizontal line 
node in the gap on a and (3 25]. In contrast, we show 
in Fig. Q] that for the non-chiral (b) pairing state the 
susceptibility for a-axis B-fields decreases approximately 
linearly towards x«/2 at T = 0, again consistent with 
the expectations of Eq. |U These results confirm that in 
spite of very complicated character of the order parame- 
ter (see ref. (25[) having zeros on the a and (3 sheets of 
the Fermi surface, the overall behavior of the calculated 
susceptibility tensor is consistent with the expectations 
from superfluid 3 He [rl H3|. 

In zero magnetic field, but in the presence of non- 
zero spin-orbit interaction A the model predicts that the 
ground state is the chiral state (a) with d-vector along 
the crystal c - axis, provided that the spin-orbit coupling 
leads to a small spin dependence of the effective pairing 
nj. Choosing U' = U n about 1% larger 
W* — is sufficient to stabilize the chi- 
ral state even for large spin orbit coupling. In contrast, 
for spin-independent interactions, U' — U, the alterna- 
tive pairing states (b) and (c) ((d) and (e) II if) are the 
ground states for any value of A < (A > 0) |26j |. 

PHASE DIAGRAM IN A MAGNETIC FIELD 

In Fig. [2] we show the free energies of the (a) (b) and (c) 
symmetry pairing states as a function of external mag- 
netic field both for JJ || c and H _L c. The model param- 
eters were chosen to make the chiral state (a) stable at 
zero field. For the field in ab-plane the free energy of the 
chiral state (a) increases with field more slowly than for 
the (b) and (c) states. This means that in large in-planc 
fields the (a) state becomes relatively more stable than 
(6) or (c). On the other hand, in an c-axis field the chiral 
phase increases its free energy much faster than the other 
phases, until at a certain critical field, B t , the (b) or (c) 
solutions become more stable. Therefore at the field B t 
we expect a "spin flop" type phase transition from a d- 
vector oriented along the c-axis to one where the d-vector 
lies in the a — b plane. For the parameter values used in 
Fig. [2] this critical field is B t « 20 mT. to one where the 
d-vector lies in the a — b plane. For the parameter values 
used in 

This prediction is consistent with what one expects 
from the analogous superfluid state in 3 He- A. In bulk 
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FIG. 2: (Color online) Condensation free energy at T = 90mK 
for three triplet order parameters of different symmetries, (a), 
(b) and (c) as a function of external field H \\ c (solid lines) and 
H _L c (dashed lines), respectively. The parameters used for 
calculations are: A = —0.02t where t = 80meV is the effective 
7 band nearest neighbor hopping integral, and U' = 1.0011I7. 



3 He there is no preferred symmetry direction in space 
and so the d-vector simply rotates continuously to re- 
main perpendicular to the applied field. But in thin films 
the d-vector is pinned, and only rotates at a finite crit- 
ical field, known as the Freedericksz transition [27j]. In 
the case of Sr2Ru04 the transition is not simply a rota- 
tion of the chiral d-vector but is also a transition from 
a chiral to non-chiral pairing state of different symme- 
try. The two distinct solutions shown in Fig. [2] (note 
that (b) and (c) are essentially degenerate) have different 
entropies and hence this Freedericksz-like spin-flop is a 
first order thermodynamic phase transition. 

One can also ask whether the specific case of Fig [2] 
is typical for a general set of model parameters. In our 
effective Hubbard model the spin-orbit interaction enters 
the Hamiltonian both directly via A in Eq. (0, and also 
in the spin dependent pairing potential U' /U. These are 
not really independent, since a full theory, such as a spin 



fluctuation model 



would include both effects 
we show a part of the full 
111 for a specific choice of 



HUME 

on the same footing. In Fig. 
phase diagram of our model 
the parameters A and U' /U. The inset shows the effect of 
the c-axis magnetic field in shifting the phase boundary 
between the chiral state (a) and the alternative (6) and 
(c) states. At any point in the phase diagram where (a) 
is stable for zero field, there is a definite critical field for 
which the spin flop to (b) or (c) takes place. However it 
is only a small region of the phase diagram, close to the 
H = phase boundary, where the transition takes place 
in fields as small as 20 — 40 mT. The implications of the 
Murakawa et al. experiment are that the parameters 
do indeed lie in this region. The alternative, is that the 
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FIG. 3: Critical interaction U' /U as a function of A in the 
presence of external field H || c at T = 0. The critical c-axis 
field of the (a) to (b) or (c) spin flop transition depends on how 
close the physical parameters lie to the zero-field (a) — (6) / (c) 
phase boundary. 
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FIG. 4: (Color online) The expected phase diagram in pres- 
ence of external field H \\ c in the region below H C 2. 



spin flop only takes place at a much larger field, and thus 
would never be observed in any c-axis field below H C 2 = 
75 mT, in contradiction to the simplest interpretation of 
Murakawa et aZ.'s result. 

If the relevant parameters do allow the d-vector rota- 
tion to occur in a relatively low field, then we expect a 
phase diagram similar to Fig. [5] This hypothesis makes 
a definite and clear testable prediction, namely the exis- 
tence of a new first order phase boundary in the H — T 
plane, which has not yet been observed. A phase dia- 
gram of the type shown in Fig. [4] is required from the 
experiments 0, 31]. The NQR experiments Q indicate 
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FIG. 5: The temperature dependence of the entropy jump 
along the transition line B — Bt = 20 mT and A = — 0.02i. 

the state with d-vector lying in the (ab) plane is real- 
ized in a c-axis magnetic field of order B = 20 mT and 
larger, and that this d-vector orientation must be present 
at these fields for all temperatures below T c . On the other 
hand the nonzero Kerr signal obtained at fields around 
10 mT [3l[ (see the points in Fig. (3a) of this reference 
and the corresponding discussion in the text) is a clear 
indication that at these fields the chiral state is stable. 
Accepting the usual interpretations of these two experi- 
ments one is forced to accept the general structure of the 
phase diagram as in Fig. 3] as the simplest picture which 
is consistent with both experiments, with a possible mar- 
gin of uncertainty of about 10 mT in the position of the 
B = B t transition line. 

An obvious question is why has the predicted phase 
transition in Fig. [4] not been observed experimentally? 
Experimental data [32| show multiple superconducting 
phases at the magnetic field parallel to the (ab) plane, 
but none have been reported for fields along the c-axis. 
The answer to this question we suggest is that the phases 
in question have essentially identical values of the quasi- 
particle energy gap, |d(k)| on the Fermi surface, and 
therefore to a very high level of accuracy they have es- 
sentially identical entropy. This is because to leading 
order only the direction of d(k) changes at each point 
on the Fermi surface, and not its magnitude. To con- 
firm this we have calculated the entropy change along 
the transition line B = 20 mT and this is shown in the 
Fig. [5J The total change in entropy is very small, of 
order 10 _6 fcs per formula unit, suggesting that it might 
not have been detected in specific heat experiments. The 
entropy change in Fig. [5] corresponds to a latent heat of 
at most 0.003mJ/mol, compared to the zero field specific 
heat jump AC/T of about 28mJ/K 2 mol at T c = 1.5K 
[22| . Assuming that lattice strain or inhomogeneity leads 
to some smearing of the ideal first order phase transition 
this small entropy change could easily be masked by the 
experimental noise. 

We should note that in the present calculations we have 
included the influence of the magnetic field in the Zeeman 
energies only and neglected the orbital effects altogether. 
This means that the vortex contribution to the condensa- 
tion energy of both the (a) and (b) phases is assumed to 
be the same. There is no reason to expect a significantly 
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FIG. 6: (Color online) Temperature dependence of the spin 
dependent interaction anisotropy ratio U' /U necessary for 
A = — 0.02t in the presence of external c-axis field (H || c) 
to stabilize the chiral phase. The symbols (a) and (6) denote 
the stability regions of corresponding phases. 
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FIG. 7: (Color online) The temperature dependence of the 
condensation free energies for states 'a' and 'b' in a c-axis 
magnetic field B = 30 mT for A = -0.02i and U' = 1.0001(7. 
The insets show the data on the expanded scale for temper- 
ature well below (upper left) and above crossing point (lower 
right). 



different vortex lattice response in the different phases, 
since the only changes are in the direction of d(k) on the 
Fermi surface and not in its magnitude |d(k)|. 

CONSTRAINTS ON MODEL PARAMETERS 

As mentioned above, the experiments 0, Hl[ indicate 
that the region where the chiral (a) state is stable is re- 
stricted to the low field part of the H — T plane be- 
low a critical field which must be in the range of about 
B t w 10 — 20 mT. This places a number of constraints 
on the possible parameter values which we can use in our 
model Hamiltonian. 

Using our model Hamiltonian, Eq. [5l we can ask what 
are the possible values of the spin-dependent pairing in- 
teraction needed to qualitatively describe both experi- 
ments. In our model this is the ratio U'/U = IP^/LP' 
for opposite spin pairing compared to equal spin pairing. 
The other model parameters are fixed beforehand, includ- 
ing the known band structure and the values of interac- 
tion parameters U\\ and U± which were previously fitted 
to other experiments [2oT |. For simplicity we have taken 
the spin-orbit coupling to have a fixed value, A = — 0.02i, 
where t — &§meV is the 7 band in-plane nearest neighbor 
hopping integral in our tight binding band structure fit. 

With all other parameters fixed, we then calculate the 
minimum spin dependent interaction enhancement U'/U 
required to stabilize the chiral state in magnetic fields of 
up to B — 20 mT. It turns out that this minimal value 
of U'/U changes with temperature, as shown in Fig. [6] 

The temperature dependence of U'/U shown in Fig. [6] 
arises because the condensation energies calculated for a 
given value of U'/U for (a) and (b) states cross as func- 
tion of temperature as shown in Fig. [7] For a small value 
of U'/U for which the (b) state will be stable at low tem- 
peratures, we find that this state is destabilized relative 
to the chiral state closer to T c . 



Let us stress that the model we are using is, in princi- 
ple at least, valid for arbitrarily large spin-orbit coupling. 
However, the parametrization of the Fermi surface which 
we have used ceases to be valid for |A| > O.K. There- 
fore in order not to refit the normal state band structure 
for every value of A we have limited the calculations to 
small values of A. In fact most of the calculations pre- 
sented in this work have been done for A = — 0.02i. This 
restricted choice does not limit the validity of our conclu- 
sions, because for every value of the spin-orbit coupling 
parameter, A, it is possible to find a corresponding value 
of U'/U which will stabilize the chiral state at T = OK 
as discussed earlier (c.f. Fig. 0]). Close enough to this 
boundary of stability a c-axis magnetic field will rotate 
the d-vector and stabilize one of the other four states 
|2rj| independently of the sign or magnitude of A. This 
minimal value of U'/U depends not only on A and B but 
also on temperature, as shown in Fig. [6l However, this 
uncertainty in the value of U'/U precludes an unique the- 
oretical determination of the phase diagram on the H — T 
plane without additional assumptions. 

Assuming a constant (i.e. temperature independent) 
value of U'/U which leads to stabilization of the chi- 
ral phase at the lowest temperatures would lead to a H- 
T phase diagram with the chiral phase occupying most 
of it. The non-chiral state would be limited to low 
temperature-high field corner. Other assumptions about 
the temperature dependence of U'/U could produce dif- 
ferent topology phase diagrams. Rather than try to ex- 
plore all these possibilities we simply assume that the 
boundary between the stable phases at low c-axis fields 
is given by a line at approximately B — B t — 20 mT 
(shown in Fig. [J) as dictated by experiment. From this 
assumption we have then calculated the corresponding 
temperature dependence of the minimal anisotropy U'/U 
which is necessary to stabilize the a-phase along the line 
(as shown in Fig. We find this approach more di- 
rect than trying to determine the expected temperature 
dependence of U'/U a priori from a spin fluctuation feed- 
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back mechanisms 33] stabilizing chiral phase. 

We note that our effective pairing Hubbard model does 
not really include the full self-consistent effects of spin 
fluctuation feedback (33|, which could make the effec- 
tive pairing interactions U and/or U' temperature de- 
pendent. However, the predicted phase transition places 
severe constraints on the superconducting pairing mech- 
anism when one takes into account a sizable spin-orbit 
coupling and the detrimental effect of the c-axis B-field 
on the chiral phase. 



SUMMARY AND CONCLUSIONS 

To summarize, we have explored the role played by 
spin-orbit coupling in determining the superconducting 
states of strontium ruthenate in the presence of c-axis 
and ab-plane oriented magnetic fields. We showed that 
the d-vector rotation can provide a consistent interpre- 
tation of both the NQR and Kerr effect data, provided 
the spin-orbit coupling constant A is small enough. Un- 
fortunately, LDA estimates of A are significantly larger 
than the values we have assumed in the calculations pre- 
sented here 



A large value of spin-orbit coupling has 
previously been cited as evidence against the d-vector 
rotation picture [3, 35, [36|. Despite these objections, 



it is clear from our results that even for large values of 
A the transition field B t for d-vector rotation, in Fig. [21 
can lie within experimental constraints 0, l3l| , but only 
if the pairing interaction spin anisotropy U' /U is such as 
to make the free energy difference between (a) and (b) 
phases finely balanced so that a small field of order 20 
mT is sufficient to cause the d-vector rotation. We have 
shown that for any value of A there is corresponding 
anisotropy U'/U such that this balance can be achieved. 
It is possible that such balancing is a consequence of a 
spin-fluctuation feedback mechanism [33|], or emerges 
directly from the full microscopic spin fluctuation theory. 
But in the absence of a detailed theory of this effect 
further progress in the field must await experimental 
measurement of the effective spin-orbit parameter for 
the quasiparticles at the Fermi surface. 
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